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Abstract 

We study theoretically the quantum dynamics of an electron in the sym- 
metric four-level double-dot structure under the influence of the monochro- 
matic resonant pulse. The probability amplitudes of the eigenstates relevant 
for the quantum dynamics are found from the solution of the non-stationary 
Schrodinger equation. The first-order correction term to the solution ob- 
tained through the rotating wave approximation is calculated. The three- 
level double-dot dynamics and the two-level single-dot dynamics, as well as 
the off-resonant excitation process, are derived from the general formulae for 
corresponding choices of the pulse and structure parameters. The results ob- 
tained may be applied to the solid-state qubit design. 

PACS number(s): 03.67.Lx, 73.21. La, 73.23.-b, 78.67.Hc 



I. INTRODUCTION 

Recently the low-dimensional semiconductor structures containing a small number 
of electrons in the size-quantized conduction band have attracted much attention. 
One of the reasons for that interest consists in their potential applicability to the 
quantum information processing [Tj. It is commonly believed that those structures 
may be scaled up to form the quantum register with appropriate number of qubits. 
During the past decade many proposals for the semiconductor qubit realization 
were made. Here we mention the phosphorous donors embedded in a silicon host 
[21 E] and a wide class of the systems based on the quantum dots (QDs) (see, e.g., 

m ei in in 13 cdi in ca ennui)- 

One of the features characterizing the coherent evolution of a quantum system is 
the Rabi oscillations induced by an external field. If the system was initially in one of 
the eigenstates of unperturbed Hamiltonian, it starts to oscillate under perturbation 
between its eigenstates giving rise to a broad class of different phenomena observed in 
the experiments. In view of the quantum algorithm realization, the Rabi oscillations 
are often considered as the process generating the desired qubit-state evolution. The 
most popular setup to study the Rabi oscillations of the electron population in the 
semiconductor nanostructures is based on the double-dot structure containing a 
single electron in the conduction band ^3 or, alternatively, on the artificial 
H2 molecule formed by the donor implantation techniques (the so-called "charge 
qubit" ^7j). In what follows we shall consider the QDs keeping in mind the general 
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character of the conclusions made below. The standard method usually exploited 
for an electron charge manipulations in such structures is based on the electrostatic 
field control. By applying an adiabatically switched voltage one can lower or raise 
the potential barrier that separates the QDs thus allowing an electron to tunnel 
between the QDs or localizing it in one of the QDs. The quantum evolution of 
an electron clearly demonstrates the Rabi-like behavior that originates from the 
coherent electron tunnelling between the QDs [TBI I16j . 

There is another scheme for the manipulations with a single electron confined in 
the double-dot structure. It is based on the resonant interaction of an electron with 
the coherent electromagnetic pulse [TH]. The laser field, instead of electrostatic one, 
is applied to the double-dot structure, and the quantum state engineering is realized 
via the optically-induced transitions between the size-quantized electron levels. As 
it was shown, the pulse parameters (the frequency, duration, and amplitude) can be 
chosen so as to drive the electron, localized initially in the ground state in one of 
the QDs, to the ground state of another QD via the state delocalized over QDs and 
used here as a "transport" state. If the states localized in different QDs are viewed 
as the Boolean states and 1, then, e. g., the electron transfer between them may 
be considered as the unitary operation NOT. The idea was initially proposed by 
Openov [IE] and then developed further in the works [IHl 1201 12U 1221 123 I2U The 
influence of strong electromagnetic fields on the tunnelling phenomena in several- 
level nanostructures was studied in Refs. [2H1 12H] • It was shown that a laser with 
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appropriate power and frequency can drive the electron between quantum wells in a 
finite-size quantum well nanostructure or localize it in one of them. In Ref. ^H] the 
opposite effect - i.e., a weak-laser-induced electron transfer between two quantum 
dots, was considered. In that model of the one-electron quantum dynamics, the 
assumption of instantaneous spreading of an excited electron over the QD structure 
was made or, equivalently, the matrix element of optical dipole transition between 
the ground and the "transport" states was supposed to be much less than the matrix 
element for tunnelling of an excited electron between the QDs. This allowed to 
describe the dynamics within the framework known in the atomic optics as the 
three-level (or A) scheme. The probability amplitudes to find an electron in the 
states localized in each of the QDs and in the "transport" state delocalized over the 
QD structure, were found as functions of the time, the pulse parameters, and the 
structure parameters. Provided that the quantum evolution is coherent, this process 
describes the three-level Rabi oscillations of the electron population. However, in 
all of these studies, the rigorous quantitative analysis of the assumption concerning 
the instantaneous spreading of an excited electron over the double-dot structure has 
not been presented. 

In this work we focus on the quantitative study of the coherent quantum dynamics 
of an electron in the symmetric four-level double-dot structure under the influence 
of the resonant laser pulse for arbitrary tunnelling rates between the excited states 
of the QDs. Here we give the detailed derivation of the analytical expressions for the 
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probability amplitudes of the electron eigenstates relevant for the quantum dynam- 
ics. The results will be presented in terms of the basis states of isolated QDs. In this 
picture, the probability amplitudes are the explicit functions of both the tunnelling 
matrix element and the matrix element of the electron optical dipole transition. We 
show that the character of the system evolution is determined by the ratio between 
these matrix elements. The three-level double-dot and the two-level single-dot dy- 
namics are derived from the general formulae as the limiting cases. We examine 
also the off-resonant excitation scheme that is very promising for the experimental 
realization of the proposed method of electron state manipulation. 

The paper is organized as follows. In Section II we present the description of the 
model and obtain the general solution for the electron dynamics in the four-level 
double-dot structure. The important particular cases of the three- and two-level 
dynamics as well as the off-resonant Raman-like excitation are considered in Section 
III. Section IV contains the results of numerical simulations. The conclusions are 
summarized in Section V. 

II. MODEL AND GENERAL SOLUTION 

Let us consider the double-dot structure (see Fig.l) containing an electron in the 
size-quantized conduction band. For the sake of simplicity, we suppose the dots 
A and B to be identical. The existence of at least two one-electron orbital states 
| A (B) 0) and \ A (B) 1) (ground and excited) in each of the QDs is assumed, with the 
one-electron wave functions <£a(b)o ( r ) — ( r I A(B)O) and <£a(b)i ( r ) = ( r I A(B)l), 
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respectively. Provided that the distance between the QDs is sufficiently large, the 
wave functions of the QD ground states are localized in corresponding QDs, and 
their overlap can be neglected. The overlap between the ground state and the excited 
state belonging to different QDs will be neglected as well: (A(B)O \ B(A) 1) ~ 0. The 
excited levels are chosen to be close to the edge of the potential barrier separating 
the QDs. They couple through the electron tunnelling [TH] . 

The Hamiltonian of an electron confined in the symmetric four-level double-dot 
structure is 

H = e (\A0) (A0\ + \B0) (B0\) + e x {\A1) {Al\ + \Bl) (Bl\) + 

+ [-V(t)\Al)(Bl\+H.c), [> 

where Eq and E\ are the one-electron energies of the ground and excited states, 
respectively (the same for both QDs); V (t) > is the matrix element for the 
electron tunnelling between the excited states of the QDs, that, in general, may be 
a time-dependent function. 

We consider the quantum evolution of an electron under the influence of the 
electromagnetic field that induces the optical transitions between the ground and 
excited states in each of the QDs (\A0) <-» \A1) and \B0) <-> It is convenient 

to examine that evolution as a complex process including both the optical excitation 
of an electron in one of the QDs and the tunnelling of the excited electron into the 
other QD. The model Hamiltonian has the form 

H = H + [A (t) (\A0) (Al\ + | SO) (Bl\) + He] , (2) 
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where A (t) is the matrix element of the electron-field interaction. In what follows we 
shall suppose the matrix elements V and A to be real and not show explicitly that 
they are time-dependent. The criterion of applicability of the model Hamiltonian 
(2) is expressed by the inequalities |A| <C u>iq, V <C u w , where u>io — S\ — e§. 
Besides, we shall suppose that there are no additional levels localized in the near 
neighborhood of the energy S\. 

The state vector of the system may be represented in terms of the eigenstates of 
isolated QDs as 

!*(*)>= E cn(t)\n). (3) 

n=A0,B0,Al,Bl 

Let an electron be initially localized in the ground state of the QD A, i.e. (0)) = 
| AO). The quantum evolution of the state vector is governed by the non- stationary 
Schrodinger equation 



.01* (*)) 



H\*(t)), 



(4) 



or, in the matrix form, 
.8 



'dt 



( c A0 \ 

cbo 
cai 

\ CBl ) 



( e A \ / c A0 \ 

e A c B0 

A E\ —V cai 

V A -V e x J V cbi ) 



(5) 



with the initial conditions c n (0) = 5„ 5 ao (hereafter h = 1). 

The straightforward diagonalization of the Hamiltonian matrix amounts to the 
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set of eigenvectors 

\l)=u-(\A0) + \B0))/V2 + v-(\Al) + \B1))/V2, 
|2) = u + (\A0) - \B0))/V2 + v+(\Al) - \B1)) / y/2, 
|3) = V .{\AQ) + \BQ))/y/2-u-{\Al) + \B1))/V2, 
|4) = w+ (|^0> - \B0))/V2-u + (\Al) - \B1))/V2, 



(6) 



and to the corresponding eigenenergies 

E 1 = e + l (co_ + yJwl+4\*) , E 2 = e + \ (u + + y^+IF) , 

£3 = eo + | (w- - y/u 2 - + *\ 2 ) , £ 4 = e + § ("+ " V<4 + 4A *) > (7) 

E 4 < E 3 < 

where m±~ — (1 — i>± ~ 1 — and u;± = c<j 10 ± V. (Here we restrict our- 

selves by the third-order terms in the parameters j^- <C 1. Note that the expansion 
of (7) over the small parameters ^-<1 gives rise to the Bloch-Siegert term ^ in 
the eigenenergies E i: i = 1 — 4.) The matrix diagonalizing the right-hand side of Eq. 
(5) has the form 



where 



Ui = ^( u - u+ ),u 2 = ±( v - v+ ) (9) 

(The columns of D are the transposed eigenvectors \k) , k — 1 — 4). 

The electron dynamics is most easily revealed in the representation of the in- 
stantaneous eigenstates, Eqs. (6), of the system Hamiltonian. In this basis, the 
time-dependent state vector of the system takes the form 

l*(*)> = EM*)l*>, ( 10 ) 
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where the instantaneous probability amplitudes (t), k = 1 — 4, are related to the 
probability amplitudes c n (t) by the matrix D: 



c(t) = D&(t) 



(11) 



where c (t) = [c A0 (t) , c B0 (t) , c A1 (t) , c B i (t)] T and a (t) = [a x (t) , a 2 (t) , a 3 (t) , a 4 (i)] J 



The Schrodinger equation in the new basis reads 
.d\® (*)> 



or, in the matrix form. 



dt 



( ai \ 
V «4 / 



D*HD-iD^\ |$(*)) 



( Ex //_ ^ 
£ 2 // H 
/il E 3 



(12) 



\ 


^ ai \ 




a 2 




a 3 


/ 


\ a 4 / 



(13) 



with the initial conditions ctx{2) (0) = u_( + )/v^2, 0.3(4) (0) = v-t+\/y/2. Here E 



E + ( 1 + ^ ) , #3(4) ~ 



and /x± 



'dt 



x 



1 - 



IX 2 



1(2) 



m agree- 



-(+)/ — (+) 

ment with the approximations for u± and v±. Besides, we have neglected the terms 

- v3 Mi ( — ) in Ei assuming them much smaller than the Bloch-Siegert term. This 



dt \ uj± 



approximation requires that 



9. (A.) 

dt \u ± J 



<C 1 and imposes the constrain on the 



speed of changes in both the tunnelling matrix element V and the matrix element 
of the electron-field interaction A. The terms ~ ^~§i (z^) m l jL± have been rejected 
as well. At this stage of consideration we see that the pairs of coefficients (ai, a 3 ) 
and (a 2 , a 4 ) evolve independently from each other. The four-level problem of Eq. 
(5) thus reduces to the two- level ones that can, in principle, be solved by one of the 
standard methods developed earlier |2Tj . 
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To proceed further we have to concretize the electron-field interaction term A in 
Eq. (2). In what follows we shall consider the situation where the tunnelling matrix 
element V is time- independent. It is worth to mention that for the adiabatically 
switched voltages the Hamiltonian matrix in Eq. (13) is diagonal and its solution 
is straightforward. As an example we give here the solution describing the electron 
behavior under the influence of the electrostatic field E 7^ E (t) = 0). The 
expression for c follows directly from Eqs. (11) and (13). The probability amplitudes 
for an electron to be in the ground states of the QDs have the form 




where A = eE J <P*a(b)o ( r ) Y{ Pa{b)i (r)dr is the matrix element of optical dipole 
transition between the states \A(B)0) and \A(B)1). The population is therefore 
localized (up to the \q/uj\q <C 1 terms) in the ground-state subspace and exhibits 
the two-level Rabi oscillations at the frequency Qq = XlV/ujf . The probabilities 
to find an electron in the states \A1) and \B1) are of order of Xl/uf <C 1 and 
oscillate at the frequency V. One sees that even for the static fields, a substantial 
electron state evolution in the double-dot structure occurs for the characteristic time 
T ~ l/^o- in principle, the electron oscillations may be utilized for the qubit-state 
engineering but this process seems to be too slow in comparison with the resonant 
optical driving JS] and unviable in view of the decoherence. However, this effect 
should be taken into account if quantum operations are performed through the 
sequence of electrostatic voltages |§J E] since in this case the electron transitions 
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between the localized and delocalized states induced by the static fields bring about 
an unwanted qubit dynamics, i.e., a computational error. 

The central part of our investigation will be devoted to the interaction of an 
electron with the time-dependent resonant pulses. For the sake of simplicity we 
consider a monochromatic square pulse of the amplitude E , the duration T, and 
the frequency uj: 

E (t) = E cos (cot) [6{t) -6{t-T)}. (15) 

In this case 

A = A cos {ut) [9 {t) -6{t-T)}, (16) 

where 6 it) is the step function. The frequency u of the laser pulse may be detuned 
from the resonant frequencies u± by the values S± = u — u±, where \6±\ C w. 

As we see from Eq. (13), it is sufficient to analyze the dynamics of just one pair 
of the coefficients, e.g., a± and a 3 (the dynamics for a 2 and a 4 is then revealed by 
the substitution u_ — > u + in the results obtained for a x and a 3 ). Transforming the 
coefficients a\ and a 3 according to the formulae a\ = aie _ ^ £o+a '-^ and a 3 = a 3 e~ i6ot 
and inserting the expression (16) for A into Eq. (13), we arrive at the set of two 
coupled linear differential equations for the coefficients d\ and a 3 (an analogous set 
is obtained for the coefficients a-i and 64): 

i~ ai = *1 cos (cot) fix + ^ (1 - e^) e" l<5 -*a 3 A' = A — (17) 

ia 3 = cos (ut) fi 3 + ^ (1 - e- i2wt ) e i5 -% ' ° °w_ ' 

Here we restrict ourselves to the terms linear on the small parameter \ /u in /i± and 
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retain the Bloch-Siegert term ~ X'^j 00 . As we shall see below, the account of this 
term is necessary for obtaining the correct result within the first-order approximation 
on the parameter \' Q /u) <C 1. 

Usually, at this point the rotating wave approximation (RWA) is made and the 
fast oscillating exponents e ±t2lJjt in Eqs. (17) are omitted. The solution thus accounts 
only of the one-photon processes conserving the energy of the system. This approx- 
imation is valid if the frequency ujiq w uj w uj± dominates the Rabi frequency that 
is of the order of |Ao|. Since this requirement on the pulse parameters is inherent to 
many optics setups, the theoretical predictions based on the RWA are in excellent 
agreement with the experimental data. However, any possible extension of the RWA 
seems to be very instructive in view of the quantitative estimate of corrections to 
the results obtained by the RWA. In Ref. [2H] the authors suggested a simple and 
clear way of how to calculate the first-order correction term to the RWA solution. 
Using the adiabatic elimination procedure for the virtual two-photon states, they 
were able to find the probability amplitudes for a two-level system beyond the RWA. 
Here we shall utilize their method to solve the set of Eqs. (17). 

According to Ref. [2HI, the coefficients a\ and 03 can be sought in the form: 

at « 4 0) + a[~ ] e- i2 ^ + a[ +) e i2 ^ , 

(18) 

a 3 « 4 0) + ajfV^' + 4 +) e i2 ^ , 
where the higher-order terms proportional to e ±t2mujt , m > 1, are dropped. From 
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the Eqs. (17) and (18) we obtain the set of six equations: 



id^ 



ia 3 
%d\ 
id 3 
%d\ 
id 3 



(°) : 

(-) 
(-) 
(+) 
(+) 



Adiabatic elimination yields 



2a;4 } + 



A'n 



(of - «i + >) 



-j(5_t„(-) 



■*(a. 



a 3 

(+) 



4°») 



,(op 



2cja 



(+) j_ A'o„iJ_t„(+) 



(19) 



,^1 ( ) r~s^ 



o, 4 +) ~ o 



(20) 



and 



4 +) 



A, 



4u- 



a< 0) e *-«, 



which in turn amount to the set 

(0) = 5A^ a (0) + Xo p -i8-t n (0) 



.(°) 



8u> "1 



5\' 2 oJO) , A', 



«3 



_tJ0) 



(21) 



(22) 



?a 3 ^ = -^fa 3 + ^e^-ai 
Note that this procedure enables us to calculate only the first-order correction term 

~ Aq/cj to the RWA since account of the higher order terms makes the set (22) 

incompatible. 

The substitution af^ = af^ e^ l ~^~ , 4°^ = d^e l ~^~ transforms the set of Eqs. 
(22) into 



;si (0) = 



(0) 



■~ (°) A'n H 



(23) 



e v /a 



2 ^ ' ' a l 

that is equivalent to the following second-order differential equation: 

~ (°) , • f x 5A 'o\ ~ (o) . A'o_( ) n 



(24) 
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with the initial conditions (0) = j^j, cii"' (0) = — «^7j- 

Of course, we have to justify the adiabatic approximation used in Eqs. (20) - (21) 
by imposing the requirement on the pulse switching time tq: 



~ (o) 



Ao 



U) 1 < T < T. 



(25) 



In what follows, however, we shall continue to handle the ramp pulses since the 
accurate calculation carried out for an adiabatically switched pulse brings about 
simple renormalization of the matrix element Ao conserving the total character of 
the ramp-pulsed dynamics (2E1- Note that the first of inequalities (25) ensures the 
applicability of Eq. (13) at t < r . Besides, we assume that r <C | Ao | - This 
unnecessary but very useful condition minimizes the influence of the pulse shape on 
the Rabi oscillation pattern. 

The solution of Eq. (24) is straightforward; transforming it back to the coefficients 



and af\ one has 



V2 



e 2 



^ cos (2f2'_t) - i^h sin (2f2'_t) 



An 



fl (0) 

"3 



V2 



e 2 



cos 



(26) 



(2ft'_*) - i S ~~^^ sin (2tt_t) 

where Q'_ = \J A'q + 5 ,2 _ 1 4, S'_ = 5_ — 5A'q/ 4c<j. It is easy to verify that these 
expressions satisfy the set of Eqs. (19). 

As a result, for the coefficients in the laboratory frame we obtain 



CAO(BO) 



ft (t) 



. 6±t 



f (t) ± / + (t) 



e 2 



cos (20tt) -i^sin (2fi±t) 



A 2 



-i2wt 



sin (2fi±t) 



CAl(Sl) 



e 

3_t 



is-t f — 



fr (t) ± e-^ft (t) 
-zsin (2Q ± t) + 
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(1 



e i2uJt ) cos (2ft±t) 



(27) 
(28) 



where 

Sl± = ^2X^ + 51 (29) 

are the Rabi frequencies and A = -^=, e± = e + u±. 

The formulae (27) - (29) describe the general type of the coherent one-electron 
evolution in the symmetric four-level double-dot structure driven by the resonant 
monochromatic pulse. One can check that the normalization condition 

E \c n {t)\ 2 = l + 0(\l/uj\5 ± /uj) (30) 

n=A0,B0,Al,Bl 

is fulfilled with the accuracy adopted through the calculations. 

We couldn't reveal any noticeable effect on the dynamics that would be produced 
by the small parameters 5±/u. They appear in the expression for the Rabi frequen- 
cies f2' ± , Eq. (26), and may compete with the terms <5±/Ag if |Ao|/w > |5±/A |. Since 
we are not interested in detailed consideration of the system dynamics in that range 
of parameters we have neglected the terms ~ 5±/u in Eqs. (27) - (29). Note that 
the terms \ 5±/uV, \q5±/uj 2 which also contain this small parameter may be com- 
parable to the terms of the order of Aq/u; 2 that have been omitted in the solution, 
and hence must be omitted as well. 

III. QUANTUM DYNAMICS IN THE STRONG AND WEAK 

TUNNELING REGIMES. 

In this Section we analyze the results obtained above for various choices of the 
pulse and structure parameters. From the general formulae (27) - (29), we derive 

15 



the expressions for the probability amplitudes corresponding to the situations where 
the characteristic tunnelling energy V is either much greater or much less than the 
matrix element of optical dipole transition |A |. Besides, we investigate also the 
electron dynamics in the strongly detuned Raman-like regime. 

A. The three-level quantum dynamics. 
If an electron being excited in one of the QDs tunnels into the other QD during a 
time much shorter than the Rabi oscillation period, one can speak about the simul- 
taneous electronic excitation in both of QDs. In other words, one of the hybridized 
states, i.e. (\A1) + \B1)) j ^2 or (|AL) - \B1)) / \pl (equally-weighted in each of the 
QDs), is excited. Both of these states are the eigenstates of the stationary Hamil- 
tonian Hq, Eq. (1). It seems then preferable to expand the state vector (t)) over 
the eigenstates of H instead than over the states of isolated QD basis. Doing so 
and using the resonant approximation we may consider only one of hybridized states 
that is formally equivalent to the setting V ^> \Xq\. The quantum dynamics of our 
system thus coincides with that of the three-level nanostructure. This situation was 
studied in the works ^H] - JS! without, however, paying enough attention to the 
mathematical proof of that proposal. In what follows we shall consider this case in 
detail and show to what extent the electron dynamics may correspond to the scheme 
just sketched. 

The condition of the strong tunnel coupling between the excited orbital states of 
the QDs, as compared to the optical dipole coupling between the ground and excited 
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orbital states of the single QD, is 

2V = \8+-8-\ > |A |. (31) 

We are interested in the resonant electron-pulse interaction, when the pulse fre- 
quency uj matches one of the resonant frequencies u>± and is strongly detuned from 
the other one. For definiteness, let the pulse frequency uj to be close to the fre- 
quency cj_ so that \5-\ <C \8 + \ and, as it follows from Eq. (31), \8 + \ 3> |A |. This 
choice of the pulse frequency corresponds to electron transition from \A(B)0) to 
(\A1) + \B1)) J \[2. Obviously, there are two different situations concerning the mu- 
tual relation between the pulse detuning 5_ and the value of |A |, i.e. <C |A | 
and \8-\ |A |. The first inequality corresponds to the resonant electron-pulse 
interaction, whereas the second one describes the off-resonant Raman-like coupling 
(see Sec. Ill C). 

Since here we consider the resonant case, all of the conditions imposed on the 
system parameters may be summarized in the following inequalities: 

|<L| < |A | < \5+\, (32) 

where 8 + — 2V and uo- = uo + — 2V. 

Making use of Eqs. (27) - (29) and taking into account that 

8+ - i6 -^( l +w\ A 2 (\8+\ \ 

cos (2Q±t) - i— =- sin (2Q±t) w e V ±/ + isgn (<5± W sin M-^-t (33) 
4sZ-t o± \ 2 I 
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at | Ao | "C \S±\, let us rewrite the Eq. (3) in the form 

|$) = I e -fcot [1 + /- (t) _ ^e" iy * sin (Vt)l | AO) + 



4V 2 

-1 + / - (t) + i^e~ iVt sin (Vt)j | -BO) + (34) 

It is easy to calculate from Eq. (34) the probability of the state inversion pbo (t) = 
\(B0 | ^ (t))| 2 after the pulse of the duration = h/Xq (the so-called 7r-pulse) is 
off: 

, . 5 2 n 2 \ 2 2 fVn\ 

PB ^ ) = l -mW-W 2 ^ [ml) (35) 

The first two terms in Eq. (35) correspond to the results of Ref. jTHj where 
the off-resonant electron transitions to the state (\A1) — \B1)) j \/2 were completely 
neglected. The third term arises due to account for such transitions. Their contri- 
bution to the electron state evolution (34) is proportional to the small parameter 
| < 1 and results in the relative phase and amplitude shifts between the coefficients 
cai and cbi that indicate on the finite tunnelling time T tunn ~ 1/V between the QDs. 
Moreover, the oscillations at the frequency 2u (so-called Bloch-Siegert oscillations) 
affect, to some extent, the ideal three-level oscillation picture. 

The results obtained show that the three-level scheme can be used for the descrip- 
tion of electron dynamics if the conditions | < 1, ^ -C 1 are satisfied. The errors 
introduced due to the presence of a nearby forth level are of the order of y. 
B. The two-level quantum dynamics. 

Next we study the opposite case of small V when a substantial spreading of an 
excited electron between the QDs occurs after many Rabi oscillations in single QD 



have completed. Such an electron dynamics is realized in the double-dot structure 
where the tunnel coupling between the excited states of the QDs is rather small as 
compared to the electron-pulse coupling: 

2V = \8+-8-\ < |A |. (36) 

The condition (36) may be satisfied for two different pulse designs, i.e. for both 
| A 1 <C \5-\ , \5 + \ and |<L.| , |<5 + | <C |A |. The first inequality corresponds to the off- 
resonant single QD excitation whereas the second one characterizes the resonant 
two- level Rabi oscillations in the same QD. The off-resonant case does not reveal 
significant two-level dynamics since an electron stays predominantly localized in the 
state \A0) (the population of the state \A1) is of the order of \q/5± <C 1). We shall 
focus our attention on the resonant transition for which the conditions 

V,|<L|,|<f+|«|A | (37) 

are satisfied. 

In this case an electron oscillates between the ground and excited states of the 
QD A for the pulse durations T <C l/V. To prove this statement we make use of 

the approximations Q + — Q- « ^f^| an< ^ e ~ %Vt ~ ^ — no ^ s f° r the ^ me 

domain t <C l/V. Inserting them into the Eqs. (27)-(29) and retaining in the time 
dependencies the terms up to the first order in Vt, we have the following expression 
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for the state vector: 



|*) = e 

+e -ie-t 



-ieot 



fo (0 + i^fe's- cos (2f2_t) 



| AO) + cos (2fi_t) \B0) \ + 



/r(i) + i^U^-^sin(2Q_0 



|A1) 



A M£ 
2^^, 2 



e 2 



sin (20_t) 

(38) 



that demonstrates the two-level electron evolution in the QD A is slightly distorted 
by the excitations in the QD B. Setting in Eq. (38) V = 0,5- = we find the 
state-vector evolution for the two-level system being in the exact resonance with the 
external pulse 



|ijr\ — e -isot 
+e -ie-t 



cos (A t/2) - i^e- i2ult sin (A t/2)] | AO) + 
-i sin (A t/2) + ± (1 - e i2 ^) cos (A t/2)1 |A1). 



(39) 



If one adopts the scheme where the quantum information is encoded into the 
ground and first excited electron (or exciton) states of the single QD ([S], [El, |2~9*]). 
the population transfer into neighboring QD should be considered as the information 
leakage from the computational subspace and the corresponding error probability 
may be evaluated with the help of Eq. (38). The computational error introduced 
by the Bloch-Siegert oscillations alone is deduced from Eq. (39). 

C. The electron excitation driven by the strongly detuned pulse. 

Finally, we shall examine the case of the Raman-like off-resonant excitation of 
an electron in the double-dot structure (|Ao| -C \8±\). This mechanism of the op- 
tical quantum-state engineering is currently under extensive investigations because 
of important properties that distinguish it from the other optical schemes (see, e.g., 
[30]). First, the excited (auxiliary) states are populated only virtually that allows 
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one to localize the electron population almost completely in the ground-state sub- 
space {\A0) , 1 50)}. Provided that the states \A0) and \B0) constitute the qubit 
basis one can therefore operate with the quantum information trapped in the logical 
subspace for any time. This, in its turn, simplifies the state evolution design and 
prevents the qubit from the decoherence induced by spontaneous photon emission 
from the excited levels. Second, the population transfer realized via the off-resonant 
excitations is quite robust against the pulse imperfections such as the uncontrollable 
detunings and the timing errors. The quantum optics provides one with wide class 
of schemes specially developed for those purposes. Quite recently several attempts 
have been made to adopt those schemes for the solid-state objects possessing of the 
atomic-like spectrum, e.g. the QDs |21| . the QDs combined with cavity QEDs jjj, 
the rf-SQUIDs HH]. 

Some features of the qubit state evolution based on one-electron quantum dy- 
namics in the symmetric double-dot structure driven by the strongly detuned pulse 
have been outlined in the work [21]. Here we consider this effect as the particular 
case of the four-level double-dot dynamics studied in Sec. II. Choosing the system 
parameters to satisfy the inequalities 

|A | < \5-\ , \8+\ < to, |A |<V, (40) 

we get from the general formulae (27)- (29) the following expression for the state 
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vector 



I*) ps e 



* £0+ 




| AO) +isin 



( 



25-5. 



V 



+ 



xh) \B0) . (41) 



The Eq. (41) describes the two-level Rabi oscillations at the frequency £1$ = 
VX 2 j25_5 + similar to those induced by the electrostatic field (see Eq. (14)). The 
Rabi frequencies of these processes are very different from each other, viz. VLq/VLs ~ 
\5-5 + \/u 2 <C 1. This makes the optically-driven oscillations more preferable for 
a qubit-state engineering due to their higher speed as compared with that of the 
electrostatic driving. As it follows from Eqs. (14) and (41) the quantum dynamics 
in both cases is frozen if V/|Ao| approaches zero. This effect may be explained in 
terms of the destructive interference between the probability amplitudes of the hy- 
bridized states represented by the symmetric and antisymmetric superpositions of 
the excited states of isolated QDs. Since those states become nearly-degenerate with 
the decrease of the parameter V/|Ao|, their probability amplitudes sum up in QD A 
whereas they cancel one another in QD B. Despite of these states are empty during 
the population transfer, the importance of their assistance to the process becomes 
more clear from this analysis. 

It is worth to note that the system evolution described by Eq. (41) cannot result in 
an arbitrary rotation of the qubit-state vector on the Bloch sphere since it contains 
only one time-dependent parameter 9 = fist corresponding to the polar angle (the 
azimuthal angle is fixed and equals to it/ 2). To overcome this obstacle one should 
break the symmetry of the structure and use at least two driving pulses with different 
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parameters to implement the desired rotation [2*U] . 



IV. NUMERICS 

To illustrate the analytical results of Sec. II and Sec. Ill we have performed 
the numerical simulations of electron dynamics in our structure. The Eq. (5) with 
A defined by the Eq. (16) was integrated within the time interval < t < 3T 
(where T = vr/Ao) for |Ao|/cu = 10~ 3 and 5_ = 0. This choice of the pulse-structure 
parameters corresponds to that usually realized in the QD systems where ujiq ~ 10 -2 
eV and |Ao| ~ 1CT 5 eV for the pulse strength Eq ~ 1 -j- 10 V/cm. Since we are 
interested in demonstration of the transition from the three-level double-dot scheme 
to the two-level single-dot scheme, the ratio V/| Ao| was varied from 0.01 to 10. 

The numerical plots showing the time dependencies of the populations p n = 
\c n \ 2 , n = AO, B0, Al, Bl are presented in the Figs. 2 (a)-(d) for V = 5; 1; 0.3; 0.05, 
respectively. For large but finite values of V/|Aq| the three-level Rabi oscillations 
picture involving the states \A0) , \B0) , (\A1) + \B1)) j \/2 becomes non-ideal due to 
the excitation of the state (\A1) — \B1)) J see Eq. (34). This effect is clearer 
seen in the representation of isolated QD basis since the phase and amplitude shifts 
between p^i and Pbi provide us with the measure characterizing the difference be- 
tween the electron populations of the excited levels in the QDs A and B. We see that 
for V/|Ao| = 5 (Fig. 2(a)) the results of Sec. Ill A may be still applied while for 
V/|Aq| = 1 (Fig. 2(b)) the regular oscillation pattern is destroyed when t > 3T and 
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the identification of excitation process is hardly possible. Figure 2(b) clearly demon- 
strates the strong modulations of the optically-induced Rabi oscillations caused by 
the electron tunnelling if V xs |A |. Further reduction of the ratio V/|A | amounts 
to the qualitative changes in the population dynamics. Figure 2(c) indicates the 
importance of the single-QD processes even for V/|Ao| = 0.3. When V/|Ao| = 0.05 
(Fig. 2(d)) we observe several almost ideal two- level Rabi oscillations in the QD A 
slightly modulated by residual dynamics in the QD B. The total reorganization of 
oscillation pattern that marks the transition from one excitation scheme to another 
occurs for V/|Ao| ~ 0.01. The reduction of the Rabi frequency by a factor of 2 
and the depopulation of the states belonging to the QD B are clearly seen from the 
numerical plots that confirms the results obtained above. 

V. CONCLUSIONS 

There are a lot of proposals for the qubit design that use the basic quantum 
properties of low- dimensional objects to encode, to process, and to store the quantum 
information. The existence of purely theoretical frameworks is of great importance 
since they allow us to capture the principal aspects of idealized evolution of the 
system under consideration and then to examine it further at the more profound 
level. An exact solution describing the qubit dynamics is often readily achieved due 
to the simplified structure of the model. It is therefore desirable to look for the 
model that would include the main features characterizing the qubit and, at the 
same time, enable the analytical treatment of the dynamical problem. 
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In this paper we have studied in detail the one-electron double-dot structure 
proposed as the candidate for a qubit implementation [TS]. The quantum operations 
in the structure may be realized by applying the resonant electromagnetic pulse 
driving an electron between the QDs. Within this scheme we have generalized the 
results recently obtained for several different pulse-structure setups [TH], [21], [2B] 
and have pointed on some delicate aspects concerning their applicability which, to 
our knowledge, had never been clarified before. As we have shown, the efficiency of 
one or another scheme is conditioned by the value of the ratio between the matrix 
element of electron tunnelling and the matrix element of optical dipole transition. 
The mathematical model of the one-electron excitation process has permitted to 
study the coherent evolution of the system beyond the rotating-wave approximation. 
The numerical results have confirmed those obtained analytically for the parameter 
choices corresponding to the three- and two-level dynamics. 

The results presented in this paper may be also applied to the two-electron sym- 
metric double-dot structure [32] and to the other systems possessing the same spec- 
tral properties, say, to the superconducting devices jSH EH|. Besides we suppose 
that the effect of the structure asymmetry on the electron dynamics [22] m &y be 
treated in the same way. 
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Fig. 1. Schematics of the states for a single electron confined in the double-dot 
structure. 
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Fig. 2. The electron populations of the ground states \A0) (solid line), \B0) (solid 
dashed line) and the excited states \A1) (thin line), \B1) (thin dashed line) of the 
symmetric double-dot structure versus the pulse duration T (in units of Aq x ) for 
the pulse-structure parameter choice <5_ = 0, |Ao|/o; = 1CT 3 and (a) V/|Ao| = 5, (b) 
V7|A | = 1, (c) V7|A | = 0.3, (d) V/\X \ = 0.05. 
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Fig. 2(b) 
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Fig. 2(c) 
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Fig. 2(d) 
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